Nonlinear Interferometry via Fock State Projection 
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We use a photon-number-resolving detector to monitor the photon-number distribution of the 
output of an interferometer, as a function of phase delay. As inputs we use coherent states with 
mean photon number up to seven. The postselection of a specific Fock (photon-number) state 
effectively induces high-order optical nonlinearities. Following a scheme by Bentley and Boyd [S.J. 
Bentley and R.W. Boyd, Optics Express 12, 5735 (2004)] we explore this effect to demonstrate 
interference patterns a factor of five smaller than the Rayleigh limit. 

PACS numbers: 42.50.St, 42.50.Ar, 85.60.Gz 



Classical optics studies light fields containing macro- 
scopic numbers of photons, while quantum optics mainly 
focuses on fields in which the mean photon number is 
much less than one. Light fields in between these two 
extremes have been relatively unexamined, because most 
detectors that are sensitive to single photons cannot dis- 
tinguish one photon from two or more. Photon num- 
ber resolving detectors have been developed, and have 
recently been applied to studies in quantum optics. In 
particular the Poisson photon number distribution of a 
coherent state has been confirmed 00] and non-classical 
photon statistics from parametric down conversion have 
been observed [3j. It should be mentioned that it is also 
possible to reconstruct the photon number distribution 
of a light field without photon number resolving detec- 
tors, either with quantum homodyne tomography 
or with avalanche photodiodes of varying efficiency 

In previous experiments, whether performed with sin- 
gle or multi-photon detectors, static field distributions 
have been measured. Here we report on the use of multi- 
photon detectors to monitor, in real time, the output 
of a Mach-Zchnder interferometer as we scan its phase. 
The inputs to the interferometer are coherent states with 
mean photon number up to seven. Projecting the output 
onto Fock states reveals a myriad of highly nonlinear re- 
sponses hidden behind the familiar, linear, measurement 
of the mean photon number. This is a canonical exam- 
ple of measurement induced optical nonlinearities, which 
have already found applications in quantum information 
processing with linear optics 00] and in a demonstration 
of a quantum controlled-NOT gate ^3 . We exploit these 
hidden nonlinearities to simulate an N-p hoton absorber, 
a substance in development 0, H2I Il3| that has gained 
interest for its application in quantum photolithography 
|l4j ] . In this manner, we implement the scheme of Bent- 
ley and Boyd |l5l ] to obtain interference patterns at one 
fifth the Rayleigh limit with classical light. 

Our detector is the visible light photon counter 
(VLPC) 0, El E3, which was developed for a parti- 
cle tracking system at Fermilab based on photon gener- 
ation in scintillating fibers. While similar in its oper- 
ation to an avalanche photodiode, it differs in two im- 
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FIG. 1: Example pulse height distribution of the VLPC. Ver- 
tical dashed lines show decision thresholds. The single photon 
count includes dark counts. 



portant respects. First, it is sensitive to photons with 
wavelengths up to ~ 30 /im, so that it must cooled to 
7K and extensively shielded from the large background 
of room temperature photons. Second, due to the de- 
tails of the VLPC design, an incident photon generates 
a tightly localized avalanche of roughly 40,000 electrons. 
Since the active area of the device is much larger than the 
size of an avalanche, photons that are absorbed in differ- 
ent locations within the detector generate independent 
avalanches. One can determine the number of absorbed 
photons in a light pulse by measuring the height of the 
current pulse from the detector. 

Figure [2 shows a representative pulse height histogram 
from the VLPC. The peaks from absorption of one, two, 
etc. photons can be distinguished. Their nonzero width 
is due predominantly to fluctuations in the number of 
electrons produced per photon. This multiplication noise 
makes the peak for absorption of n photons y/n broader 
than the single photon peak, because it is the result of n 
independent processes. 
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FIG. 2: Experimental Setup. The waveplate, birefringent 
crystal (biref), and polarizer form a polarization Mach- 
Zehnder interferometer. The output state is collected into 
a single mode fiber and sent to the VLPC. 



In the multi-photon resolving experiments mentioned 
above, the photon number distribution was determined 
by sending all current pulses to a multichannel analyzer 
that measured the pulse height histogram. The photon 
number distribution was then determined by fitting the 
histogram to a sum of Gaussians. This is successful, but 
the use of a multichannel analyzer limits the allowable 
data rate. Also, the determination of how many photons 
were absorbed is made only in the fitting stage, after 
the experiment is over. We employed a different method 
that allowed us to determine, on a shot by shot basis, 
how many photons were detected. We used several com- 
parators (Analog Devices AD96687) that compare the 
amplified VLPC current pulse to seven variable thresh- 
old levels Lfc. We first measured a pulse height histogram 
by scanning two closely spaced levels and monitoring the 
rate of pulses between them. We then set the thresholds 
of the comparators to the minima between the photon 
peaks, as shown by the vertical lines in Fig. ^ A pulse 
that lay between levels and Lk+i was recorded as the 
detection of fc photons. Due to fluctuations in the number 
of electrons per avalanche, it is possible that fc photons 
can be recorded as k + 1 or k — 1. This becomes more 
of a problem at higher k, as the peaks overlap more. In 
the experiments reported here this overlap is reasonably 
small for k < 5. If more accuracy is required, this can be 
corrected if one has all the nonzero count rates . 

Our optical setup is shown in Figure El Our 
light source was a 780 nm wavelength, 100 fs pulsed 
ThSapphire laser. Because our detector begins to sat- 
urate at around 5 MHz, we reduced the repetition rate 
R lcp of the laser to 40 kHz with a Q-switched amplifier 
system. The laser light was heavily attenuated with a 
neutral density filter and sent into a polarization Mach- 
Zchndcr interferometer. The phase delay of the inter- 
ferometer was changed by tilting a birefringent crystal 
inside it. The light from one output was then collected 
into a single mode fiber and sent to the VLPC. 

For our purposes [lj| , we can take the laser light to be 
in a coherent state \a), given by [2fJ 
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where \k) is a fc-photon Fock state and a is the dimen- 
sionless electric field amplitude of the coherent state. The 
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FIG. 3: (Color online) (a) through (g): Count rates vs. phase 
for the detection of one to seven photons from an input state 
with no = 3.95. (h): average photon number computed from 
weighted sum of count rates. Solid lines show theoretically 
expected values. 



mean photon number of such a state is (n) — \a\ 2 . The 
interferometer transforms the input state \a) into the 
output state |cvsin (0/2)), so that the intensity of the 
output varies sinusoidally with phase. Losses in the sub- 
sequent optics and detector transform the interferometer 
output into the detected state \r) ■ a sin (0/2)), where r] 2 
is the total detection efficiency of the setup. Thus the 
detected state is equivalent to the output of an interfer- 
ometer fed by a weaker input state. 

Figures|2Ia) to (g) show the count rates for a maximum 
average detected photon number n max of 3.95. Each plot 
shows the count rate for detection of k photons, or equiv- 
alently the overlap of the detected coherent state with the 
fc th Fock state. An important feature of each count rate 
is the non-sinusoidal behavior. For photon numbers k 
less than n max there is a local minimum at the position 
of maximum light intensity. Ideally, the fc-photon count 
rate Rk is given by 



R k = i? rep |(fc|asin(0/2))| 5 



-rw sin 2 (0/2) (rWxSin 2 (0/2)) A 
fc! 



(2) 
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FIG. 4: Single photon count rate vs. interferometer phase at 
several mean photon numbers: (a) 0.463, (b) 0.99, (c) 1.99, 
(d) 3.95, (e) 7.18. 
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FIG. 5: Calculated phase sensitivity vs. phase when measur- 
ing the average photon number, and the overlap of the state 
with the one, two, or four photon Fock states for n max = 4. 



As the intensity of the light increases, its overlap with the 
fc-photon state increases only until (n) = k. Because Rk 
is maximal when (n) = k, the number of count rates that 
peak before 4> = 7r is an indication of n max . To determine 
this value quantitatively, we fit the single photon count 
rate to equation[21 which yielded n max = 3.95. Using this 
value we generated the lines for R2 to R7 and (n) . Note 
that i?7 is significantly higher than predicted from equa- 
tion[2J because eight photon absorption is non-negligible. 

As can be seen from their complicated phase depen- 
dence, the individual count rates are composed of many 
frequency harmonics. However, when they are combined 
as X)fe=o kpk as in Figure Qh, only the fundamental fre- 
quency remains, ft is this sinusoidal sum that is usually 
associated with the output of an interferometer. By mon- 
itoring the variation of the photon number distribution, 
the rich underlying behavior is brought to light. 

To observe the appearance of a dip in the single pho- 
ton count rate, we took five scans with different values of 
"max- Figure 0] shows the observed single photon count 
rates. At low photon number, the single photon rate 
follows the sinusoidal intensity pattern, but as n max in- 
creases, it smoothly deforms into a double peak structure. 

Given that i?7, for example, is much more sharply 
peaked than (n) , and therefore contains many higher fre- 
quency harmonics, one might expect that an interferom- 
eter can be made more sensitive by monitoring i?7 rather 
than (n). This is not the case. To prove this, we need 
a result from quantum estimation theory |2l| . When 
measuring a quantity O to estimate a parameter </>, the 
uncertainty in the estimate is quantified by 
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where A is the standard deviation in the measurement of 
a quantity. The standard method of interferometry is to 
measure n. For a coherent state with mean photon num- 
ber n max , A(j> = l/y/n max cos(</3/2). The best sensitivity 
is obtained at <j> = 0, where the intensity is zero. 

The analysis for measuring a fc-photon count rate is 
similar. The number of counts we get is given by a bi- 
nomial distribution, and the variance in the count rate 
is (ARk) 2 = Pfe(l — pk)/N, where N is the number of 
laser pulses used, and pk is the probability that the state 
projects to |fc). This leads to 
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The phase setting that gives optimum sensitivity is now 
a function of n max and fc. Sensibly, the worst sensitiv- 
ity occurs when (n) = n max sin 2 (0/2) = k, where Rk is 
maximal and thus insensitive to changes in the output 
intensity. 

Figure |3] compares the calculated sensitivity of these 
two measurement schemes. The usual intensity measure- 
ment is significantly more sensitive than monitoring a 
fc-photon count rate, except for the single photon count 
rate, which yields equal sensitivity at <f» = 0. Here the 
output intensity is nearly zero, there are never any pulses 
with two or more photons, and the single photon rate is 
proportional to the intensity. 

The presence of these higher harmonics is the result 
of the nonlinearity of the detection process, i.e. the 
multi-photon absorption probability is a highly nonlin- 
ear function of the light intensity. Materials that ab- 
sorb nonlinearly could prove useful for many tasks such 
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FIG. 6: Seven photon count rate vs. interferometer phase, 
(a) alone, and superimposed with copies shifted by (b) n, (c) 
27r/3, and (d) 2n/5. Note the change in scale in part d. 



by 2tt/N is a measure of the harmonics at and above N. 

In conclusion, we have presented real-time measure- 
ments of a photon number distribution using a detector 
that can resolve photon numbers up to seven. We used 
this technique to measure the output of an interferom- 
eter, which showed that the photon count rates behave 
highly nonlincarly and vary on a sub-wavelength scale. 
While the phase sensitivity of an interferometer cannot 
be improved by this method, we use it to demonstrate 
sub-wavelength interference patterns using only classical 
states of light and the inherent nonlinearity of detection. 
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as photolithography. A recent proposal by Boto et. al. 
|l4j suggests a way to obtain interference features a fac- 
tor of N smaller than the Rayleigh limit. It requires 
both an iV-photon absorbing resist and light in the en- 
tangled state |A)|0) +e jAr *|0)|A), where the kets denote 
the two interfering light modes, and <j> — 2irx/\ is the 
phase difference between them at a position x on the re- 
sist. In this case the probability of absorbing N photons 
goes as cos 2 (A</>/2), which yields interference fringes N 
times finer than the Rayleigh limit. Later, Bentley and 
Boyd dl| demonstrated that entanglement is not essen- 
tial. Their scheme obtains interference fringes a factor 
of N smaller than the Rayleigh limit using an M photon 
absorber and classical light. The major drawback is that 
for high visibility fringes, M must be much greater than 
N. The scheme was presented as follows: interfere two 
laser beams to produce a standard interference pattern 
on the photoresist, and increase the relative phase of the 
two beams by 2tt/N for each pulse. The M-photon ab- 
sorption is proportional to J2 i= i(EiE*) M , where Ei is 
the electric field on the resist from pulse i. Shifting the 
phase of the pulses averages out the spatial frequencies 
on the resist below the desired frequency of A/TV. 

Here we present a similar but experimentally easier to 
realize scheme. Figure HJa) shows the 7-photon absorp- 
tion vs. phase at n roax = 3.95. Panels (b), (c), and (d) 
show the results of a data post-processing in which the 
same pattern is superimposed with copies shifted by mul- 
tiples of 2"7r/2, 2"7t/3, and 2ir/5, respectively. The essence 
of the proposal now becomes clear: the multiphoton ab- 
sorption peaks on the resist are narrow, and more than 
one of them will fit in the space of one wavelength. The 
variation in the count rate that remains under shifting 
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